Duality Between Spatial and Angular Shift in Optical Reflection 
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We report a unified representation of the spatial and angular Goos-Hanchen and Imbert-Fedorov 
shifts that occur when a light beam reflects from a plane interface. We thus reveal the dual nature of 
spatial and angular shifts in optical beam reflection. In the Goos-Hanchen case we show theoretically 
and experimentally that this unification naturally arises in the context of refiection from a lossy 
surface (e.g., a metal). 
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Introduction. — In the 17th century Newton was the 
first to surmise that the center of a reflected beam should 
present a small spatial shift A in the plane of incidence, 
relative to its geometrical optics position fl*]. More than 
two centuries afterwards, in 1947 Goos and Hanchen 
(GH) 2| were able to quantitatively measure such a shift 
(see [3| for a literature survey since that time). The GH 
shift is typically in the sub-wavelength domain; it has 
become technologically important in recent years since it 
directly affects the modes of optical waveguides and mi- 
crocavities [4, 5] and has great potential for (bio)sensor 
applications [6] . Theoretically, the GH shift has been ex- 
plained at various levels and several generalizations have 
been discovered [7, 8,, 9} . Amongst the latter it was pre- 
dicted that the axis of the reflected beam should display 
a small angular deviation O from the law of specular re- 
flection 6>inc = ^ref [I, IQ] • Interestingly, it took about 50 
years since the original experiment performed by Goos 
and Hanchen to actually observe such angular shift in 
the microwave [11] and the optical [3] regimes. Presently, 
it is common wisdom that spatial and angular GH shifts 
are two different phenomena observables in two mutually 
exclusive regimes: the spatial GH shift occurs in total 
reflection (reflected intensity = incident intensity) [l2| . 
while the angular GH shift appears in partial reflection 
(reflected intensity < incident intensity) [10]. 

In this Letter we show that this separation is artiflcial. 
We present a unifled description for the spatial and angu- 
lar GH shifts that will appear as two aspects of a unique 
beam-propagation phenomenon. We show that such du- 
ality between spatial and angular shift is rather general 
and also applies to the Imbert-Fedorov (IF) effect, which 
is a shift normal to the plane of incidence [13] that has 
drawn considerable interest lately 0, 0, &]. Finally, 
for the GH shift we show that uniflcation naturally arises 
in the context of reflection from lossy surfaces. For this 
case we also furnish an experimental demonstration that 
the spatial and angular shifts occur simultaneously. 

Our Letter is structured as follows: We give flrst a 
qualitative picture of the envisaged uniflcation. Then, 
we furnish a rigorous theoretical analysis of the beam- 
propagation problem and show that the unifled descrip- 



tion actually holds for both the GH and the IF shifts. 
Finally, we demonstrate, both theoretically and experi- 
mentally, that for the GH case lossy reflecting surfaces 
naturally induce uniflcation. 

Qualitative picture. — Consider a system consisting of 
two homogeneous, isotropic media of dielectric constants 
ei and £2, filling the half-spaces z < and 2; > 0, re- 
spectively, as shown in Fig. 1. A monochromatic beam 
of light of wavelength Aq and waist wq propagates along 
the central wave vector ko in the region z < before im- 
pinging upon the plane interface of equation z = that 
separates medium 1 from medium 2. Detailed deriva- 
tions of the angular and the spatial displacements O and 
A for this system have already been reported elsewhere 
0, B @, El, E, [195 and will not be repeated here. We 
merely quote the basic results in the form 



Ax = Xolm[Dx], 6> 



9l/2)Re[D, 



(1) 



where Aq = Ao/(27r), and 60 = 2Xo/wo is the angular 
spread of the incident beam [20] . The expressions above 
are valid for both the GH and the IF shifts where the 




FIG. 1: (color online). Scheme of the beam reflection at the 
plane interface. Here ko and ko are the central wave vectors 
of the reflected beam as predicted by geometrical and wave 
optics, respectively. The waist of the incident beam is located 
at the origin of the laboratory Cartesian frame xyz. 



2 



coefficient Dx is equal to 
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in the GH case, and to Dx = 2i[(rp + rs)/rx]cot6 in 
the IF case |21|. Here = Re[Dx], D\ = lm[Dx], 
and the index A is a label for the two linear polarizations 
parallel (A = P or TM) and perpendicular (A = 6' or TE) 
to the central plane of incidence x-z. Moreover, rx = 
rx (0) = i?A exp(z(/)A) is the Fresnel reflection coefficient 
[22| evaluated at the central angle of incidence 6>, where 
Rx = \rx\ and (/)a = argrA- 

From Fig. [1] it follows, using elementary geometri- 
cal considerations, that the total beam displacement ob- 
served at distance / from the origin is expressible as a 
linear combination of Aa and 6a (supposedly 6a <C 1) 
of the form: 



Sx{l) = Ax + iex. 



(3) 



From now on, in order to avoid trivial repetitions, we 
will consider explicitly the GH case only. However, all 
our conclusions for the GH case can be straightforwardly 
extended to the IF case with minimum effort. 

Equation ([3|) supports the hypothesis that Aa and 6a 
are two different manifestations of a unique phenomenon, 
as they can be connected by such an elementary geomet- 
ric relation. From Eq. ([2]) it follows that in the case of 
total reflection {Rx = 1), the reflection coefficients re- 
duce to a pure phase factor rx = exp(z^A), thus = 
and the reflected beam undergoes a purely spatial shift 
Aa- Vice versa, if the reflection coefficients are strictly 
real, as in the case of air-to-glass partial reflection, then 
(j)x G {0, ±7r} ^ D\ = 0, and the beam undergoes a 
purely angular deflection 6a. However, when the re- 
flection coefficients are complex^ then both Aa 7^ and 
6a 7^ 0, and a unified description of spatial and angular 
GH shifts becomes mandatory. Since for a plane inter- 
face between two lossless media the Fresnel coefficients 
are either purely real or pure phase factors [22[ , then we 
have either Sx = I 6a or ^a = ^a, respectively, and spa- 
tial and angular GH shifts are mutually exclusive. On 
the other hand, when either one or both of the media 
are lossy, the separation between spatial and angular GH 
shifts becomes artificial, and unification naturally occurs. 

Theoretical description. — As a ffist step towards uni- 
fication, in this section we aim to deduce Eq. (|3]) from 
a perfectly general ab initio calculation, without restric- 
tions on the media forming the interface or on the trans- 
verse shape of the light beam. 

It is well known that under minimal hypotheses [1, Q 
it is possible to obtain a valid angular spectrum repre- 
sentation for the electric field of the reflected beam of 
the form A(X, Z) = (27r)-i // A(k, Z) e''^'^ d^n, where 
A(k, Z) is uniquely determined by the angular spectrum 



of the incident beam and the Fresnel reflection coeffi- 
cients [l9|. A Cartesian reference frame attached to the 
reflected beam with a scaled coordinate system is uti- 
lized in which X = /cqx^, Y = /cq^^, Z = /cq^^, and 
/co = |ko|, as shown in Fig. ([1]). Moreover, we have 
defined X = Xx"^ + and = k - • k) = 

Kix'^ + n2y'^ is the transverse part of the unit wave vec- 
tor k = k//co = + n^z^ with respect to = ko/Zco, 
with = (1 — • A^)-^/^. Here ko = ko — 2z{z • ko) is 
the central wave vector of the reflected beam as ruled by 
geometrical optics. The angular spectrum in the plane Z 
is determined by its value at Z = via the relation [iot : 



A(a^, Z) = A(a^, 0) exp(-iHZ), 



(4) 



where H = —1^3 is the so-called optical Hamiltonian that 
governs the well-known Hamilton equations of motion for 
light rays in vacuum [2^: 
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The analogy between Eq. (|4]) and the expression for the 
time evolution of the wave function of a quantum sys- 
tem in the Schrodinger picture, suggests the use of an 
enlightening quantum-like notation [3, 0, by writ- 
ing A(a^,Z) = (a^|A(Z)) and A(X, Z) = (X|A(Z)), 
where \k) = \k,i^K2) and |X) = \X^Y) are the ba- 
sis vectors in the transverse momentum and position 
space, respectively. In Eq. (|4]) the longitudinal coor- 
dinate Z has the role of a dimensionless time, then we 
can write \A{Z)) = exp(— z7YZ)| A(0)), where H is the 
Hamiltonian operator defined via (a^|7Y|a^') = —(1 — k ■ 
K — k'). In quantum mechanics unitary evolution 
implies that probabilities are conserved along with prop- 
agation, namely (A(Z)|A(Z)) = (A(0)|A(0)). In our 
case, this means that the flux of the electric field energy 
density across any plane Z = const., is independent from 
Z, namely JJ \A{X, Z)\^ dXdV = // | A(a^, 0) |^ d^,^ = 
1, where we have renormalized the field amplitude of 
the reflected beam to ensure (A(0)|A(0)) = 1. At 
any given coordinate Z the electric field energy density 
|A(X, Z)| gives the spatial beam proffie in the obser- 
vation plane X-Y. The Z-dependent centroid of such 
energy distribution (X) (Z) = // X| A(X, Z) l^dAdF, 
measures the deviation of the beam axis with respect 
to the central axis z^ 0, defined by geometri- 
cal optics. If we define the transverse position and 
momentum operators such that X|X') = X'|X') and 
K|a^') = a^'|a^'), respectively, then the centroid of the 
beam can be evaluated as (X) (Z) = (A(Z)|X| A(Z)) = 
(A(0)|e^^^Xe-^^^|A(0)) = (A(0)|X//(Z)| A(0)), where 
1^h{Z) is the position operator in the Heisenberg picture: 



(a^|X^(Z)|a^O 
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The first term on the right side of Eq. (|6]) coincides with 
the momentum-representation of the Z-inde pen dent po- 
sition operator X in the Schrodinger picture [27| . There- 
fore, from the very definition of (X) (Z), it follows that 
the expectation value (A(0)|X| A(0)) = (X) (0) gives 
both the GH and IF spatial shifts {Xh{0)) and {Yh{0)), 
respectively. The second term, which is linear in Z, is 
proportional to right side of the second equation in ([5j) 
which determines the direction of propagation of classi- 
cal rays of light, since k/k^s = tan'^^i^r^ + tan'i^2y^- As 
we expect small angular deviations from geometrical op- 
tics predictions, we can write tani^^ ^ i^i {i = 1,2), 
and the exact relation between wave and geometrical 
optics established by Eqs. dMS]) allows us to identify 
d{XH{Z))/dZ ^ ^1 with the angular GH shift, and 
d{YH{Z))/dZ ^ 1^2 with the angular IF shift. Such iden- 
tification becomes even more clear by noting that from 
[K,7Y] = and Eq. 1^ it follows that 

(X)(Z) = <X) + Z(K(1-K-K)-^/^), (7) 

which reduces to (X) (Z) ^ (X) + Z(K) for < 1 , 
and the angular brackets indicate expectation values with 
respect to the state |A(0)). Equation ([7j) has the same 



form (X) (Z) = A + Z0 as Eq. (|3|), where A depends 
on the position operator X, and on the momentum 
operator K solely, thus defining unambiguously both a 
spatial and an angular vector shift of the beam equal to 
A = (X) (0) and S = d (X) {Z)/dZ, respectively. Note 
that the dependence of Eq. ([7j) on the characteristics 
of the refiecting surface and on the polarization of the 
incident beam is contained in the form of the state | A(0)). 

Equation ([7j) establishes the first part of the announced 
unification by reproducing Eq. ([3]) that was naively de- 
duced on the ground of simple geometric reasoning. We 
emphasize that, because of its vector form, Eq. ([7j) de- 
scribes both GH and IF shifts. The next step is to demon- 
strate that refiection from lossy surfaces induces simul- 
taneously both spatial and angular GH shifts. 



Loss-induced unification. — Consider again the system 
shown in Fig. 1. Assuming air as medium 1, we can 
write the dielectric constant of medium 2 as £2 = Sr-\-iSi. 
Then, we can distinguish that medium 2 is a dielectric, 
namely £^ > 1; or a metal with Sr < 0. For both cases, 
from Eq. ([2]) and the well-known expressions for the Fres- 
nel refiection coefficients [22j, it follows that 



Dp 



2 sin 61 



+ (1 - £r) + iSi (1 - 2£r) 



[sr - sin^ e) + (ef - e2) cos2 6* + ie^ (1 - 2e^ cos2 9) ' 
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(8) 



To elucidate the role of the losses, we expand Eq. (|8|) 
in a Taylor series around = 0, and keep terms up to 
the first order in £i. For sake of simplicity, from now on 
we will consider only the metal case for which the Taylor 
expansion furnishes 

sin l9 (1 - 2 sin"^ - Sr (sin^ (9 + cos^ 0) 



Up —tri 



(sin^ e 
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(sin^ — Sr) -\- s'^ cos^ 



-2sr sinO 



sm 
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sm 



Sr COS' 



(9) 
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for P polarization, and = 6^2 sin ^/ (sin^ — €r 

D^g = —2 sinO/ \/ sin^ — for S polarization. From 
the expressions above we can see that for an ideal loss- 
less metal {si = 0) we have that Dxjs purely imaginary 
and only the spatial shift occurs [29]. However, when 
losses are "turned on" by letting < <C 1, a first- 
order real part must be added to causing the si- 
multaneous existence of both spatial and angular shifts. 
Such loss-induced coexistence between A and B is clearly 
illustrated in Fig. 2 where Eq. (j3|) is displayed for reflec- 
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FIG. 2: (color online). 6p{1) - 6s{l) at / = 11.9 cm for a 
fundamental Gaussian beam with wavelength Ao = 820 nm, 
waist Wo = 32 /xm, incident on a gold surface with S2 — 
-29.02+i 2.03. Blue line = Ap-As; green line = /(Gp-Gs); 
red line = 5p — 5s] black triangles = measured 5p —5s- 



tion at an air-gold interface. The red curve represents 
the total beam shift at distance / from the origin and it 
is given, according to Eq. (j3|), by the sum of the green 
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FIG. 3: (color online). Layout of the experimental set-up. 



curve (angular shift), and the blue one (spatial shift). 
This completes the second part of our unification pro- 
gram. 

Experiment. — The experimental set-up is sketched in 
Fig. 3. The 820-nm output of a super-luminescent light- 
emitting diode (SLED) is spatially filtered by a single- 
mode optical fiber (SMF) to select a TEMqo mode. This 
is then collimated by a microscope objective and sent 
through a Glan polarizing prism (PBS) to fix its polar- 
ization to P. A polarization modulator (PM) switches 
the beam polarization between P and S to exploit the 
fact that both spatial and angular GH shifts are polar- 
ization dependent. Then, a lens is used to focus the beam 
to the desired spot size in front of the mirror. While the 
spatial GH shift is unaffected by such beam focusing, the 
angular GH shift depends on the beam angular aperture. 
The reason for this behavior is evident from the Oq fac- 
tor in Eq. ([T]), and its physical origin is explained in 
[3]. Finally, the difference-signal from a quadrant detec- 
tor (QD) is fed into a lock-in amplifier in order to detect 
the beam displacement in the plane of incidence when 
polarization is switched between S and P. See [H, Is^l for 
further details. 

The results of the measurements are shown in Fig. [2] 
along with theoretical prediction. The good agreement 
confirms the simultaneous occurrence of the spatial and 
the angular GH shifts in the lossy regime. 

Conclusions. — In this Letter we have presented a uni- 
fied description for spatial and angular Goos-Hanchen 
shifts occurring in light beam reflection from lossy sur- 
faces. Such description applies to the Imbert-Fedorov 
effect as well. The unification theory has been worked 
out for a weakly- absorbing metal and the corresponding 
GH shift has been observed experimentally in reflection 
from an air-gold plane interface. 
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